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18 Abstract. We present observations of the infrared radiative cooling by carbon dioxide (CO 2 ) and 

19 nitric oxide (NO) in Earth’s thermosphere. These data have been taken over a period of 7 years 

20 by the SABER instrument on the NASA TIMED satellite and are the dominant radiative cooling 

21 mechanisms for the thennosphere. From the SABER observations we derive vertical profiles of 

22 radiative cooling rates (W m' 3 ), radiative fluxes (W nf 2 ), and radiated power (W). In the period 

23 from January 2002 through January 2009 we observe a large decrease in the cooling rates, 

24 fluxes, and power consistent with the declining phase of solar cycle 23. The power radiated by 

25 NO during 2008 when the Sun exhibited few sunspots was nearly one order of magnitude 

26 smaller than the peak power observed shortly after the mission began. Substantial short-term 

27 variability in the infrared emissions is also observed throughout the entire mission duration. 

28 Radiative cooling rates and radiative fluxes from NO exhibit fundamentally different latitude 

29 dependence than do those from CO 2 , with the NO fluxes and cooling rates being largest at high 

30 latitudes and polar regions. The cooling rates are shown to be derived relatively independent of 

3 1 the collisional and radiative processes that drive the departure from local thermodynamic 

32 equilibrium (LTE) in the CO 2 15 pm and the NO 5.3 pm vibration-rotation bands. The observed 

33 NO and CO 2 cooling rates have been compiled into a separate dataset and represent a climate 

34 data record that is available for use in assessments of radiative cooling in upper atmosphere 

35 general circulation models. 

36 

37 1. Introduction 

38 The terrestrial thermosphere and mesosphere have been the least explored regions of 

39 Earth’s atmosphere. They are too high for in-situ measurements from aircraft or balloons. Rocket 

40 or ground-based measurements on a global basis are not practical and do not provide the suite of 

41 measurements required for a complete characterization of these regions. In the middle 1990s 

42 remote sensing technology became sufficiently advanced to enable comprehensive satellite 
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1 observations of these regions of the atmosphere. The NASA Thermosphere -Ionosphere- 

2 Mesosphere Energetics and Dynamics (TIMED) mission was developed to explore the Earth’s 

3 atmosphere above 60 km altitude and was launched in December 2001. A fundamental goal of 

4 the mission is to quantify the energy budget of mesosphere and thermosphere. One of four 

5 instruments on the TIMED mission, the Sounding of the Atmosphere using Broadband Emission 

6 Radiometry (SABER) instrument, was specifically designed to measure the energy budget of the 

7 mesosphere and lower thermosphere [Mlynczak, 1996; 1997]. 

8 The SABER instrument is a 10-channel limb scanning radiometer. It scans the Earth’s 

9 limb continuously recording profiles of infrared radiance (W m' 2 sr' 1 ) from the atmosphere in 

10 discrete spectral intervals [ Russell et al., 1999]. The specific wavelength bands observed by 

1 1 SABER were chosen so that a variety of data products could be retrieved or derived, including 

12 kinetic temperature, ozone, water vapor, carbon dioxide, atomic oxygen, atomic hydrogen, rates 

13 of energy deposition, rates of energy loss, and rates of radiative heating and cooling. To this end, 

14 SABER measures two of the three key infrared emissions that govern radiative cooling of the 

15 atmosphere above 100 km: nitric oxide (NO) at 5.3 pm [Kockarts, 1980] and carbon dioxide 

16 (CO 2 ) at 15 pm [Curtis and Goody, 1956]. The third emission, the fine structure line of atomic 

17 oxygen (O P) at 63 pin [Bates, 1951] is not measured by SABER but has been observed by the 

18 CRISTA instrument that flew twice in the 1990s on the Space Shuttle [Grossman et al., 2000] 

19 and by numerous rocket experiments [Grossman and Vollmann, 1997] and even by high altitude 

20 balloons [Mlynczak et al., 2004], The SABER dataset is the first global, long-term, and 

2 1 continuous record of the NO and CO 2 emissions from the thermosphere. 

22 The SABER dataset is now seven years in length and has already provided some basic 

23 insight into the heat budget of the thermosphere on a variety of timescales. Mlynczak et al., 
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1 [2003; 2005; 2007a] have reported the observed nature of the nitric oxide emission in the 

2 thennosphere in response to geomagnetic storm events, demonstrating a “thermostat” effect in 

3 which the infrared emission rapidly and efficiently removes storm energy deposited into the 

4 thennosphere and thus returns it to its original state in a short period of time. The long-term 

5 effects of the declining phase of solar cycle 23 were clearly evident in the NO emission in less 

6 than 5 years [Mlynczak et al., 2007b] and 9-day periodicities in both the NO and CO 2 emissions 

7 have been related to the occurrence of coronal holes on the Sun [. Mlynczak et al., 2008], 

8 The purpose of this paper is to further examine variability in the NO 5.3 |im and CO 2 15 

9 |im emissions with the SABER dataset. In particular, the extent of solar minimum conditions 

10 during the end of solar cycle 23 offers a unique opportunity to study the radiative cooling in the 

1 1 thennosphere under exceptionally quiescent conditions. We also examine the infrared emissions 

12 as a function of latitude and altitude. In addition, the newest version of the SABER data set 

13 (vl.07) is used in this study. Lastly, the results presented here form a climate data record for the 

14 thermosphere useful for assessing general circulation models of the upper atmosphere. The 

15 cooling rate data may be obtained from the authors upon request. 

16 In the next section we review the data processing methodology. In Section 3 we present 

17 the results in the order of the processing: vertical profiles of radiative cooling (W m' ), fluxes of 

18 infrared emission (W m' ), and radiated power (W). Discussions on the fundamentally different 

19 roles of cooling by NO and CO 2 are given. In Section 4 we discuss the sensitivity of the results to 

20 parameters in the SABER retrieval algorithms and in Section 5 we give a summary. 

2 1 2. Derivation of Infrared Parameters 

22 The fundamental SABER measurement is the infrared radiance (W m' sr" ) measured in 

23 the limb viewing geometry as a function of tangent altitude. The instrument scans the 


3 



1 atmosphere from approximately 400 km tangent altitude down to 5 km below the hard Earth 

2 surface, recording an infrared radiance sample every 0.4 km. The NO 5.3 pm emission is 

3 observed as high as 300 km during periods of intense geomagnetic storms. The CO 2 15 pin 

4 emission is typically observed up to 139 km. Approximately 1400 profiles of infrared emission 

5 are obtained for each of the 10 SABER channels in one day. During the limb scan between 100 

6 and 200 km (the nominal range reported in this paper) the TIMED spacecraft moves 

7 approximately 1 degree of latitude. The nominal instantaneous field of view of the instrument is 

8 2 km at the tangent point on the Earth’s limb. 

9 The TIMED satellite is in an orbit inclined 74 degrees with respect to the equator. The 

10 SABER instrument views perpendicular to the spacecraft velocity vector. In order to keep the 

1 1 SABER instrument within its allowable range of operating temperatures, the TIMED spacecraft 

12 executes a rotation about its yaw axis every 60 days to keep the Sun from illuminating the 

13 instrument’s thermal radiator. While in the “northward” viewing mode the SABER instrument 

14 views from approximately 83°N to 55°S. After the yaw maneuver the instrument views from 

15 approximately 55°N to 83°S. This operational scenario results in the regions between 55°N and 

16 55°S being viewed continuously while the regions above 55 degrees latitude are viewed 

17 continuously for 2 months every other yaw cycle, i.e., two months of viewing followed by a 2 

18 month gap, followed by 2 months of viewing, etc. This process repeats exactly the same 

19 temporal gaps in both hemispheres every year. 

20 In the results below for zonal mean cooling rates and fluxes we present annual averages. 

21 Because of the SABER viewing geometry, the averages are over 12 months only between +/- 55 

22 degrees in latitude (corresponding to 82% of atmospheric area). In the regions poleward of 55 

23 deg we average the six months of available data every year. While not truly annual, these 6- 
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1 month averages cover locally two months during winter, two months during fall, one month 

2 during spring, and one month during summer. There is no obvious discontinuity in the cooling 

3 rates and fluxes shown in the figures below and we believe the results are substantially 

4 representative of the annual behavior, although comparisons with models are required to verify 

5 this assertion. More than 50,000 individual cooling rate profiles are observed in a 6-month period 

6 poleward of 55 degrees. 

7 The derivation of the infrared radiative cooling rates (W m' 3 ) is accomplished by 

8 different techniques for the NO and CO 2 emissions owing to the different opacities of these 

9 molecules in the limb view, and to account for the different physical processes that result in 

10 emission and cooling of the atmosphere by these molecules. To further study the energy balance 

11 we vertically integrate each profile of radiative cooling to obtain the flux (W m' ) of energy 

12 exiting the thermosphere. The fluxes are then zonally integrated to provide an estimate of the 

13 equator-to-pole distribution of radiative cooling which is central to the large-scale atmospheric 

14 circulation. The fluxes are then integrated with respect to area around a latitude circle, and then, 

15 from pole-to-pole, to estimate the total power (W) radiated each day by the thermosphere. These 

16 analysis procedures with SABER data are given in Mlynczak et al., [2005, 2007a]. We will 

17 review these here since there are fundamental differences in the approaches for analyzing the 

18 (more) optically thick 15-pm band of CO 2 and the optically thin 5.3 |im band of NO. 

19 2,1 NO(u) at 5.3 |im 

20 Infrared emission from NO is a substantial cooling mechanism throughout the 

21 thermosphere. Due to the low atmospheric density, the vibration-rotation bands of the NO 

22 molecule depart from local thermodynamic equilibrium (LTE) [e.g., Funke and Lopez -Puertas, 

23 2000]. The NO vibrational states are excited primarily by collisions with atomic oxygen with a 
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1 rate coefficient of 4.2 x 10" 11 exp (-2700/T) cm 3 s' 1 [Hwang et al., 2003], where T is the kinetic 

2 temperature. By comparison, excitation by the absorption of terrestrial radiation is several orders 

3 of magnitude smaller. Nitric oxide is also fonned through a series of chemical reactions that 

4 excite high-lying vibration-rotation bands especially during geomagnetic storms. These have 

5 been observed near 195 km by the CIRRIS-1A instrument that flew on the Space Shuttle [Dothe 

6 et al., 2002; Sharrna et al., 1996]. 

7 The approach to deriving the NO cooling rates from SABER observations has been 

8 discussed by Mlynczak et al. [2005], and is reviewed here. The emission from nitric oxide is in 

9 the weak line limit of radiative transfer. The equation of radiative transfer describing a transition 

10 in the weak line limit is given in the limb view by the expression: 

11 R(H a ) = h£X.^ V(x)dx (1) 

4 7T 

X 

2 1 

12 where R(H 0 ) is the measured spectrally integrated radiance (W m' sr' ) at tangent altitude H 0 , h 

13 is Planck’s constant, c the speed of light, v the frequency of the transition (wavenumbers, cm' 1 ), 

14 V(x) is the volume emission rate of photons at the tangent altitude, and dx is a path element 

15 along the limb line of sight. An Abel or geometric inversion may be applied to the measured 

16 radiance profile, assuming spherical symmetry about the tangent point, yielding a vertical profile 

17 of the rate of emission of energy per unit volume (W m' ). For optically thin transitions such as 

18 those of NO, the energy escapes the thermosphere completely: half of the energy is emitted to 

19 space, half to the atmosphere below. 

20 As the SABER instrument is a filter radiometer, the initial result from the Abel inversion 

21 process is a cooling/emission rate that is weighted by the filter bandpass of the instrument. To 

22 obtain the total rate of emission, a correction factor is applied to the measured in-band emission 

23 rate to obtain the emission rate for the entire band. The correction factor varies substantially with 


6 



1 


altitude for NO and has been computed from extensive model calculations discussed in Mlynczak 

2 et al. [2005]. The correction factor was verified through observations of NO spectra made by the 

3 MIPAS instrument [ Gardner et al., 2006] on the EnviSat satellite. Surprisingly, Mlynczak et al. 

4 [2005] found relatively small variation in this correction (“unfilter”) factor from quiescent to 

5 geomagnetically disturbed conditions. 

6 In the results presented below we treat the entire emission from NO as cooling the 

7 atmosphere. Terrestrial (“earthshine”) emission is insignificant as an excitation mechanism of 

8 nitric oxide vibrations, relative to collisions with atomic oxygen. The results from Dothe et al. 

9 [2002] indicate that emission from high-lying NO vibrational levels is still well below that of the 

10 fundamental (t) = 1 to t) = 0) band. The correction factor applied to the SABER measured in- 

1 1 band emission rate accounts for out-of-band contributions of the 2-1 and 3-2 vibrational states as 

12 well. Potential refinements to this process will be discussed in Section 5. 

13 2.2 C0 2 (v 2 ) at 15 pm 

14 In this section we outline the method to compute the cooling rates for the CO? 15 pm v 2 

15 fundamental band transition, with upper state designation CO 2 (01 1 0). Later in the paper we will 

16 examine the role of cooling by the first hot band in this mode, i.e., the C0 2 (02 0 — » 01 0) 

17 transition. Above 100 km, the C0 2 15 pm fundamental band is not exactly optically thin, but 

18 neither is it very optically thick. Because of this, derivation of the cooling rates for the CO? 

19 molecule is not as straightforward as for NO. Specifically, the absence of true weak line radiative 

20 transfer precludes application of an Abel or geometric inversion to simply derive a 

21 cooling/emission rate profile. In addition, radiative exchange with other layers of the atmosphere 

22 occurs, although absorption of radiation from the lower atmosphere is relatively small [Lopez- 

23 Puertas and Taylor, 2001]. Only emission resulting from collisional excitation results in a 
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cooling of the atmosphere. In addition, due to the opacity, not all emission escapes to space or 
the atmosphere below, unlike NO. We must therefore employ a technique to compute the true 
cooling of the atmosphere given these constraints. 

The SABER experiment measures emission in the V 2 bending mode of CO 2 at 1 5 pm for 
the purpose of retrieving the kinetic temperature of the atmosphere. The 15 pin bands are also 
responsible for the radiative cooling of the atmosphere. As with NO, the vibration-rotation bands 
of carbon dioxide depart from local thermodynamic equilibrium (LTE) in the upper mesosphere 
and thermosphere. The SABER temperature derivation includes detailed modeling of the 
collisional and radiative processes that drive the observed transitions from LTE [Mertens et al., 
2001]. Specifically, the Curtis matrix approach [ Lopez-Puertas et al., 1986a, b] is employed in 
the SABER non-LTE modeling and temperature derivation process and directly yields the 
infrared radiative cooling rates in Kelvin per day. These cooling rates computed at each iteration 
step are output and stored as standard data products upon convergence of the temperature 
retrieval. 

To derive the local rate dQ/dt of radiative cooling in W m" 3 for CO 2 at a given altitude we 
apply the first law of thermodynamics: 


dQ 

dt 



( 2 ) 


C p is the heat capacity at constant pressure, p the density, and dT/dt the derived radiative cooling 

rate in Kelvin per day. From the ideal gas law this can be written as: 

dQ = P_^dT_ (3) 

dt T R dt 


In Equation 3, R is the gas constant. This approach to computing cooling rates and fluxes was 
successfully applied in the stratosphere by Mlynczak et al., [1998.] In our present analyses the 
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1 pressure, temperature, and dT/dt are obtained from the SABER vl.07 dataset. Although the 

2 effect is small, we also include the variation of the ratio C p /R with altitude to be consistent with 

3 that used in the SABER data processing, thus accounting for the effects of atomic oxygen on C p . 

4 3. Results 

5 3.1 Cooling rates for NO at 5.3 |im 

6 We begin by examining the cooling rates (W m' 3 ) for NO over the seven years of data to 

7 date. The SABER instrument began routine operations in mid-January 2001. To maintain 

8 SABER’s orientation on the “cold” side of the spacecraft, the TIMED spacecraft undergoes a 

9 1 80-degree yaw maneuver every 60 days, as dictated by the progression through local time in the 

10 satellite’s 74 degree inclined orbit. Thus in a year the spacecraft executes a series of six yaw 

1 1 maneuvers, allowing the same progression of local time to be repeated each year. In all figures 

12 presented below, unless noted otherwise, we display annual averages for each of the seven “yaw 

13 years” that have been observed to date. Each “yaw year” runs from mid-January to the next mid- 

14 January. So the year 2002 runs from January 2002 to January 2003, and so on. 

15 Shown in Figure 1 are the zonally averaged, annual mean cooling rates in W nf for NO 

16 observed by SABER from January 2002 through January 2009, i.e., yaw years 2002 through 

17 2008. There are two features immediately obvious from the data. First, the maximum cooling 

18 rates occur at high latitudes near the poles, implying a strong geomagnetic influence. Second, the 

19 magnitude of the cooling rates visibly decreases each year in the series from 2002 through 2008. 

20 This decrease is concurrent with the declining phase of solar cycle 23. At the altitude of the 

21 cooling rate peak poleward of 60 degrees, the cooling rates decrease by a factor of -4.5. The 

22 cooling rate is observed to decrease throughout the entire range from 100 to 200 km over the 7 

23 years of observations. By 2008 there is very little emission observed by SABER above 180 km. 
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For carbon dioxide we find a different picture of cooling rates in the plots shown in 

2 Figure 2. Whereas the NO cooling rates showed a clear equator-to-pole enhancement throughout 

3 the thermosphere, CO 2 exhibits a relatively small gradient in equator to pole cooling above 105 

4 km. Below 105 km there is evidence of strong cooling near the poles and also at the equator. The 

5 equatorial enhancement in CO 2 cooling is visible in all 7 years of data. Similar to the NO 

6 cooling, a decrease in the strength of the CO 2 cooling over the 7 years is also visible in Figure 2. 

7 Figures 1 and 2 demonstrate a decrease in radiative cooling by NO and CO 2 in the 

8 thennosphere over the last seven years. This decrease is consistent with a decrease in the 

9 temperature of the thermosphere which affects both NO and CO 2 , although the temperature 

10 sensitivity of the NO emission is substantially larger than for CO 2 . The strength of radiative 

11 cooling is proportional to exp (AE/keT) where AE is the photon energy, ke is Boltzmann’s 

12 constant, and T is the kinetic temperature. Since AE for NO is about 2.8 times larger than for 

13 CO 2 , NO emission has a substantially larger sensitivity to a specific change in temperature. 

14 Decreases in atomic oxygen over the solar cycle also contribute to the decline in infrared cooling 

15 from both molecules. 

16 The NO cooling will also be affected directly by a decrease in the NO density that would 

17 be expected to occur during the declining phase of the solar cycle, while the CO 2 concentration 

18 has only a small anticipated source due to continued anthropogenic buildup and virtually no 

19 known sinks that are due to the solar cycle. Thus the larger observed changes in NO cooling as 

20 compared with CO 2 are due to a combination of a larger sensitivity of the emission to 

21 temperature and to an overall decrease in the NO abundance. The influence of solar cycle on 

22 radiative cooling is best shown in Figures 3a and 3b which depict the annual, global average 

23 cooling in 2002 through 2008, for NO and C0 2 , respectively. For NO in Figure 3a we have 
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1 extended the plot up to 250 km altitude, illustrating the extent of nitric oxide emission during 

2 solar maximum (and also the sensitivity of the SABER instrument.) The NO cooling shows a 

3 continued and marked decrease for each year in the sequence. The largest decrease is in the 

4 global average cooling by nitric oxide, approximately a factor of ~6.5 at the peak. In contrast, the 

5 CO 2 cooling has decreased by about 35% over the solar cycle. The smaller sensitivity of the CO 2 

6 cooling to the solar cycle is due to the fact that the C0 2 emissions originate in the lower (more 

7 dense) atmosphere, to the smaller sensitivity of the C0 2 emission to temperature changes, and to 

8 a lack of si nk s of C0 2 abundance. 

9 3.2 Fluxes of Exiting Longwave Radiation 

10 The next step in the process of assessing the radiative cooling is to integrate vertically the 

1 1 cooling rate profiles and obtain the flux of longwave radiation that exits the thermosphere. We 

12 choose to call this the ‘exiting’ longwave radiation (ELR) in analogy to the ‘outgoing’ longwave 

13 radiation (OLR) used in studies of the Earth’s tropospheric climate. The difference is that the 

14 OLR all leaves the planet whereas only about half of the ELR leaves the planet while the rest is 

15 absorbed in the atmosphere below the thermosphere. The range of vertical integration is 100 to 

16 200 km for NO and 100 to 139 km for C0 2 . 

17 Shown in Figures 4 and 5 are polar stereographic plots of ELR for NO in the northern 

18 and southern hemispheres, respectively. The scale on the color bar runs from 0.0 to 1.0 

19 milliwatts per square meter. These ELR plots show the annual average fluxes of radiation exiting 

20 the thermosphere between 2002 and 2008 and correspond to the same time periods as the cooling 

21 rates shown in Figure 1. Clearly visible are larger fluxes near the poles in all years, with fluxes in 

22 2002 and 2003 peaking in excess of 0.8 milliwatts per square meter. Also evident is a dramatic 

23 decrease in flux from 2002 to 2008. In 2002 in most of each hemisphere the ELR exceeds 0.3 


11 



1 milliwatts per square meter. In 2008 the ELR over most of the globe is well below 0.07 

2 milliwatts per square meter. Thus the effect of the declining phase of solar cycle 23 is clearly 

3 evident in the ELR. 

4 Figures 6 and 7 show the ELR for the northern and southern hemispheres from 2002 to 

5 2008 for CO?. The scale runs from 0.8 to 2.4 milliwatts per square meter. The CO 2 fluxes are 

6 substantially larger than those for NO and are reflective of the role CO? plays in cooling the 

7 lower (and more dense) thermosphere. Evident in these figures is a general decrease from 2002 

8 to 2008. However, as evidenced by the emission rates of energy in Figures 2 and 3b, the decrease 

9 in the ELR for CO 2 is much less than that for NO, on a percentage basis. As with the NO 

10 emission, the CO 2 ELR is always largest in polar regions, reflecting the high latitude 

1 1 enhancements in CO 2 cooling shown in Figure 2 in the lower thermosphere. The tropical 

12 enhancement in cooling by CO 2 shown in Figure 2 is also evident in the fluxes in Figures 6 and 

13 7. The tropical enhancement in fact exhibits 4 alternating maxima and minima in each year, 

14 perhaps implying a dynamical (tidal?) effect. The tropical maximum decreases substantially from 

15 solar maximum in 2002 to solar minimum in 2008. The extent of large fluxes of ELR from CO 2 

16 also decreases at the pole, but there is still substantial cooling, although over a reduced area, in 

17 2008. 

18 3.2.1 Zonal Mean Fluxes of Exiting Longwave Radiation 

19 Next we zonally average the ELR for NO and CO 2 shown in Figures 4 through 7 to 

20 display the annual average ELR as a function of latitude. In Figures 8 and 9 we show the annual 

21 average ELR for each year 2002 through 2008 for NO and CO 2 , respectively. The ELR is plotted 

22 versus the sine of latitude so that the abscissa is proportional to atmospheric area. The first 

23 feature that stands out from examination of both figures is that starting in 2003, the ELR 
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1 decreases each year compared with the previous year. The decrease is generally uniform with 

2 latitude although the ELR in the high latitude northern hemisphere in 2003 is comparable to that 

3 in 2002, and again comparable in 2004 and 2005. The ELR for NO decreases by about a factor of 

4 8 in the low to mid-latitudes, while the ELR for CO 2 decreases by about 30% over the same 

5 range of latitudes. The general decrease in ELR from 2002 to 2008 is again an indication of the 

6 effects of the declining phase of solar cycle 23. However, the fact that the ELR decrease is not 

7 identical at all latitudes indicates that the effects of solar variability are not necessarily uniform 

8 throughout the thermosphere. 

9 A second feature that stands out upon examination of Figures 8 and 9, and which is 

10 hinted at in Figures 1 and 2, is that the equator-to-pole gradient of the ELR is different for NO 

1 1 and CO 2 , with the gradient for NO being larger than for CO 2 . The equator-to-pole gradient of net 

12 heating (solar heating minus radiative cooling) of the atmosphere is a major factor that drives the 

13 large-scale circulation of the thermosphere, and the atmosphere in general. The fact that the 

14 gradient is larger in NO and becomes lesser over the solar cycle as both temperature and the NO 

15 abundance decrease implies a potential link between chemistry, dynamics, and radiation 

16 governed by the abundance of NO. 

17 Shown in Figure 10 is the zonally averaged ELR for the sum of NO and CO 2 . As shown 

18 for the individual ELR terms, the combined ELR is larger in all years at high latitudes than at 

19 low latitudes in both hemispheres. Because the ELR is the vertically integrated radiative cooling 

20 rate, the results illustrated in Figure 9 demonstrate that the high latitude and polar thermosphere 

21 is cooled more strongly by infrared radiation than at low and tropical latitudes. This result is 

22 fundamental. Larger radiative cooling at high latitudes than at low latitudes is known to exist in 

23 the mesosphere [. Mlynczak , 2000] and the stratosphere [Mlynczak et al., 1998], However, in the 


13 



1 troposphere, the opposite is true, i.e., the outgoing longwave radiation is larger at the equator 

2 than at the poles [ VonderHaar and Suomi, 1969]. The results of this and previous studies show 

3 that the upper atmosphere cools in a fundamentally different manner than the troposphere. 

4 Another property of the ELR evident from Figures 10 is that in addition to the general 

5 decrease in ELR over the 7 years of data, the equator-to-pole gradient of radiative cooling has 

6 weakened substantially over this time period. While the decrease in ELR implies that the 

7 thermosphere has cooled, and in the case of NO, that the abundance of NO has also likely 

8 decreased, the weakening of the equator-to-pole gradient of cooling implies that the large-scale 

9 dynamics may have also weakened considerably during this time. These infrared emissions thus 

10 provide solid evidence as to changes in the thennal structure, photochemistry, and dynamics of 

1 1 the thennosphere, and hint strongly at a coupling between NO chemistry and the large-scale 

12 dynamics and transport in the thermosphere. 

13 3.3 Global Radiative Power 

14 The last parameter that we will calculate is the global radiative power (Watts) emitted by 

15 C0 2 and NO. These have been previously been shown by Mlynczak et al. [2007b, 2008] in 

16 examining long and short term variations in the energy balance of the thermosphere. We obtain 

17 the total global power by zonally integrating the fluxes of ELR shown above, and then 

18 integrating the power in each latitude bin from pole-to-pole. To compute the power poleward of 

19 55 degrees in the hemisphere not being observed, we follow Mlynczak et al. [2005] in assuming 

20 the ratio of the power between the equator and 55 degrees latitude and between 55 degrees 

21 latitude and the pole is the same in both hemispheres. This process provides a measure of the 

22 total global power radiated by NO and C0 2 on a daily basis. The power for the full 7 years of 

23 data is shown in Figure 11. What is evident is the overall decline of radiated power consistent 
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1 with the declining phase of solar cycle 23. In addition there is substantial short-term variability in 

2 both the NO and CO 2 emissions, that has previously been linked to geomagnetic phenomena 

3 [Mlynczak et a/., 2008]. Large increases in radiative power associated with extreme geomagnetic 

4 events (e.g., the October 2003 Halloween superstorm event) are evident. In addition, in 2008, the 

5 values of NO power reach very low values, nearly a factor of 10 lower (when averaged over a 

6 60-day yaw cycle) than at the beginning of the mission. During 2008 the Sun exhibited no 

7 sunspots for over 200 days. We suggest that the SABER data over this time period offer an 

8 excellent resource for studying the influence of the Sun on the climate of the upper atmosphere. 

9 The CO 2 cooling rates, fluxes of ELR, and radiative power shown above are all computed 

10 for the fundamental band of the CO 2 molecule. In the SABER temperature retrieval process more 

1 1 than 20 vibration-rotation bands (including isotopic bands) are considered and the non-LTE 

12 problem is solved for each band, including cooling rates in Kelvin per day. To verify that the 

13 fundamental band dominates the CO 2 cooling we also computed cooling rates (W m" ), fluxes of 

14 ELR, and radiated power (W) for the first hot band of the C0 2 bending mode. Shown in Figure 

15 12 is a time series of the power radiated by the CO 2 fundamental and first hot bands. It is clear 

16 from this figure that the fundamental band dominates the CO 2 15 pm thermospheric cooling and 

17 that the first hot band is about 3.5% of the fundamental band emission. However, the decrease in 

18 NO emission over the solar cycle is such that, between 100 and 200 km, the radiated power from 

19 the first hot band of CO 2 rivals that of NO in 2008, on an annual average basis, although the 

20 emissions do peak in substantially different regions of the atmosphere. Table 1 lists the annual 

21 average power for CO 2 (fundamental and first hot bands) and NO for the 2002 through 2008 

22 “yaw years.” 

23 4.0 Sensitivity to non-LTE parameters in the SABER algorithms 
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3 

1 A key question in the analysis and derivation of the CO 2 cooling rates in W m' is the 

2 extent to which the results are sensitive to the non-LTE processes and the modeling of these 

3 processes in the SABER temperature, pressure, and cooling rate (K/day) derivation in the 

4 operational SABER data processing algorithms. We demonstrate that the derived cooling rates, 

5 in W m' 3 , as a function of altitude, are essentially insensitive to the non-LTE parameters in the 

6 SABER algorithms, because the rate of emission is essentially constrained by the SABER 

7 radiance measurements. Furthermore, the ELR fluxes and radiated power are similarly 

8 insensitive since they are derived directly from the cooling rate in W m' 3 . 

9 The primary process by which CO 2 cools the atmosphere, analogous to NO, is by 

10 radiative emission subsequent to collisional excitation by an oxygen atom. In 1970 P. Crutzen 

1 1 postulated that this process would be important in cooling the thermosphere, and since that time, 

12 substantial efforts have been expended both in the laboratory [e.g. Castle et al., 2006] and in 

13 analyzing atmospheric observations [e.g., Rodgers et al., 1992; Lopez-Puertas et al., 1992] to 

14 determine the rate coefficient for the process 

15 O + C0 2 (d) -> CO 2 + O + 667.5 cm' 1 (4) 

16 Literature values for this rate “O/CO 2 ” coefficient range over a factor of 4 from approximately 6 

17 x 10" 12 cnrV \Sharma and Winter steiner, 1990] to ~ 1.4 x 10‘ 12 cm’s" 1 [ Khvorostovskaya et al., 

18 2002]. The SABER temperature algorithms incorporate the Sharma and Wintersteiner value for 

19 this rate coefficient. 

20 In the SABER algorithms, the modeling of the non-LTE processes depends on the 

21 provision of the O/CO 2 rate coefficient, the atomic oxygen density, and the CO 2 abundance. 

22 Atomic oxygen is provided from the NRL-MSIS model [Picone et al., 2002] and the CO 2 

23 abundance is provided from the Whole Atmosphere Community Climate Model (WACCM) 
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1 model [R. Garcia, National Center for Atmospheric Research, private communication, 2006]. 

2 With these provisions, temperature is retrieved and constrained by the measured radiance at 15 

3 pm by the SABER instrument [ Mertens et al., 2001]. 

4 To assess the effect the O/CO 2 rate coefficient on the derived cooling rates in W nf from 

5 CO 2 , we ran 7 complete days of the vl.07 operational SABER algorithm in which we reduced 

6 the rate coefficient by a factor of 4 uniformly at all altitudes. This is also equivalent, in the non- 

7 LTE modeling, to varying the atomic oxygen concentration by a factor of 4, or the product of 

8 atomic oxygen and the rate coefficient by a factor of 4. We computed the cooling rates (W in’ ), 

9 ELR fluxes (W m' ) and power (W) for each of the test days and compared with the original 

10 values. The results showed that the cooling rates (W in' ), fluxes, and power did not show much 

1 1 variation despite a huge variation in the key non-LTE rate coefficient. There are large changes in 

12 temperature, pressure, density, and cooling rates in Kelvin per day, but not in the cooling rate in 

13 W m' . An example is shown in Figure 13. 

14 The explanation for the relative insensitivity of the cooling rate (W m' ) is that decreasing 

15 the O/CO 2 rate coefficient results in a decrease in the population of excited CO 2 molecules in the 

16 algorithm. However, the SABER algorithm must compensate to match the measured radiance. It 

17 does this primarily by increasing the temperature and pressure. An increase in temperature 

18 increases the population of vibrationally excited CO 2 molecules by the Boltzmann factor as a 

19 consequence of detailed balance. The algorithm thus works to provide the correct number of 

20 excited CO 2 molecules to create the emission rate necessary to match the measured radiance. 

2 1 Thus, the cooling from CO 2 in W nf is essentially unchanged. Our approach is found to be quite 

22 robust and in fact provides the cooling rates (in W m' ) essentially independent of the non-LTE 

23 model parameters used to produce the SABER data products. The previously mentioned fact that 


17 



1 


radiative excitation from the lower atmosphere is not very large also contributes to this result. 

2 We note the NO cooling rates in W m' are also derived independent of any modeled non-LTE 

3 processes by virtue of the weak line inversion and because the NO emission is essentially 

4 generated solely by collisions with atomic oxygen. 

3 

5 The relative insensitivity to non-LTE processes of the derived cooling rates in W m" 

6 presented here is a substantial result. The accuracy of the derived cooling rates is then tied 

7 directly to the absolute calibration of the SABER instrument through the measured radiances. 

8 The cooling rates also represent a data set that can be used to test directly cooling rate 

9 parameterizations in upper atmosphere general circulation models [e.g., Qian et al., 2009]. The 

10 model parameterizations are very sensitive to the non-LTE processes including the rate 

1 1 coefficient for energy transfer between atomic oxygen and carbon dioxide because they have no 

12 constraint (e.g., a measured radiance). Through such comparisons the long-standing 

13 discrepancies between renderings of the key non-LTE parameters may be resolved. The relative 

14 insensitivity to the large change in the rate coefficient also suggests that the dataset of cooling 

15 rates can form a climate data record [ National Research Council, 2004] for the thermosphere. 

16 5. Discussion and Summary 

17 We have presented 7 years of observations of the radiative cooling in the thermosphere as 

18 observed by the SABER instrument on the TIMED mission. The cooling rates in W m" are 

19 derived for emission from CO 2 at 15 pm and NO at 5.3 pm. The emissions exhibit both short and 

20 long-term variations that illustrate the sensitivity of the thermosphere on time scales ranging 

21 from daily to the 1 1 year solar cycle. The cooling by NO is substantially more variable than that 

22 of CO 2 due to a greater sensitivity to temperature and to a much larger variability of NO density 
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1 with solar activity. Decreases in atomic oxygen between 2002 and 2008 likely contribute to the 

2 observed decrease in both NO and CO 2 cooling. 

3 The variability of the cooling has several fundamental consequences. First, the larger 

4 equator-to-pole gradient in cooling by NO, and its variability, implies a potential link to 

5 dynamics and transport as the equator-to-pole gradient in net heating drives the large-scale 

6 circulation. The observed weakening of the equator-to-pole gradient in radiative cooling over 

7 these seven years is a strong indication that the large-scale thennospheric dynamics have also 

8 weakened. In addition, the larger cooling at the poles than at low latitudes in the thennosphere is 

9 consistent with the observed cooling in the mesosphere and stratosphere, but is opposite of what 

10 is observed in the troposphere. Thus the entire upper atmosphere behaves in a fundamentally 

1 1 different manner than the troposphere with regards to radiative cooling as a function of latitude. 

12 We have also shown that the cooling rates derived herein are essentially insensitive to the 

13 parameters and rate coefficients used to compute the non-LTE populations of the carbon dioxide 

14 molecule. This is because the retrieval algorithm must match the measured radiance, and thus 

15 will produce enough excited CO 2 molecules to accomplish this, whether by increasing 

16 temperature or pressure (or both) to accomplish the result. Further study on the sensitivity to the 

17 carbon dioxide abundance used in the retrieval is warranted. We would expect perhaps only 

18 minor sensitivity to CO 2 abundance given the relative insensitivity to a factor of 4 change in the 

19 primary process responsible for cooling, and that we do not expect the CO 2 used in the SABER 

20 retrieval to be off by a factor of 4. At this time we estimate the cooling rates for CO 2 emission to 

21 be accurate to better than 15%. 

22 In the case of NO, the assumption is that the observed emission is essentially all cooling, 

23 and that the Abel inversion and correction factors applied to determine the total band emission 
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1 are sufficient. There are uncertainties in the derived cooling, most notably in the process used to 

2 derive the full band cooling from the measured in-band emission. This correction (“unfilter”) 

3 factor is essentially a function of the rotational temperature of the NO vibration-rotation bands. 

4 For the vl.07 SABER data we do not make a correction for the rotational temperature as a 

5 function of time. As the thermosphere cools, we would expect the rotational temperature to 

6 decrease and thus there is a possibility that the correction factor to generate the total cooling rates 

7 would also change. If that were the case, the correction factor would in general decrease, 

8 implying that the changes with solar cycle may be larger than presented above. These 

9 possibilities will be investigated in future studies by the SABER science team. We also estimate 

10 the cooling rates due to NO to be accurate to better than 15%. This estimate is based on the 

11 accuracy of the SABER radiometric calibration (1%), the possible effects of rotational 

12 temperature uncertainties (~ 10%), and on emissions that may not be properly accounted for in 

13 correcting for the spectral response filter function of the SABER NO channel (~ 10%). 

14 Assuming these uncertainties are uncorrelated, the root-sum-square of these is conservatively ~ 

15 15%. 

16 In closing we suggest that the dataset of radiative cooling of the thermosphere by NO and 

17 COt constitutes a first climate data record for the thennosphere. The length of the data record 

18 and the apparent lack of dependence of the data products on model parameters results in an 

19 accurate set of parameters that can be used to conduct fundamental tests of general circulation 

20 models of the thermosphere. The cooling rates, radiative fluxes, and radiated power can all be 

21 compared against model calculations. This should enable tests of model chemistry and physics 

22 and in principle, resolve long-standing issues with regards to the parameterization of radiative 

23 cooling in these models. As the TIMED mission continues, these data derived from SABER will 
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1 become important in assessing long term changes due to the increase of carbon dioxide in the 

2 atmosphere. 
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2 

3 

4 Table 1. Annual mean, global power (W) radiated from the thermosphere (100 to 200 km) by the 

5 carbon dioxide fundamental band, the carbon dioxide first hot band, and by nitric oxide, from 

6 2002 through 2008. 

7 


Year 

CO 2 Fundamental 
Power (10 11 W) 

CO 2 First Hot Band 
Power (10 11 W) 

NO (At) = 1) 
Power (10 11 W) 

2002 

9.19 

0.334 

2.26 

2003 

8.85 

0.332 

1.94 

2004 

8.05 

0.306 

1.21 

2005 

7.69 

0.306 

0.947 

2006 

7.46 

0.288 

0.583 

2007 

6.97 

0.267 

0.407 

2008 

6.90 

0.220 

0.337 
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1 Figure Captions 

2 

3 Figure 1. Zonal average, annual mean cooling rates (W m' ) for nitric oxide from 2002 through 

4 2008. 

5 

6 Figure 2. Zonal average, annual mean cooling rates (W m' 3 ) for carbon dioxide from 2002 

7 through 2008. 

8 

9 Figure 3a. Global annual average cooling rate profiles, solar maximum (2002) and solar 

10 minimum (2008) and all years in between, for NO, derived from the cooling rates in Figure 1. 

11 

12 Figure 3b. Global annual average cooling rate profiles, solar maximum (2002) and solar 

13 minimum (2008), for CCF derived from the cooling rates in Figure 2. 

14 

15 Figure 4. Annual average exiting longwave radiative flux (mW m' ) for nitric oxide, northern 

16 hemisphere, 2002 through 2008. 

17 

18 Figure 5. Annual average exiting longwave radiative flux (mW m-2) for nitric oxide, southern 

19 hemisphere, from 2002 to 2008. 

20 

21 Figure 6. Annual average exiting longwave radiative flux (mW m' ) for carbon dioxide, northern 

22 hemisphere, from 2002 to 2008. 

23 
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2 

1 Figure 7. Annual average exiting longwave radiative flux (mW m' ) for carbon dioxide, southern 

2 hemisphere, 2002 to 2008. 

3 

4 Figure 8. Annual mean, zonal average exiting longwave radiative flux (W m' ) for nitric oxide, 

5 2002 through 2008. 

6 

7 Figure 9. Annual average, zonal mean exiting longwave radiative flux (W m' ) from carbon 

8 dioxide at 15 pm, 2002 through 2008. 

9 

10 Figure 10. Annual average, zonal mean exiting longwave radiative flux, (W in'"), for nitric oxide 

1 1 and carbon dioxide, 2002 through 2008. 

12 

13 Figure 11. Daily global power (W) emitted from the thennosphere by CO 2 (red) and NO (green) 

14 between 100 km and 200 km altitude. 

15 

16 Figure 12. Daily global power (W) computed for CO 2 fundamental and first hot bands in the 

17 thennosphere between 100 km and 140 km. 

18 

19 Figure 13. Results of SABER retrieved data products (radiative cooling rate in K/day; radiative 

20 cooling rate in W m' ; kinetic temperature; pressure; and density), from the nominal retrieval 

2 1 algorithm (green) and from a sensitivity test in which the quenching rate of vibrationally excited 

22 carbon dioxide by atomic oxygen was reduced uniformly by a factor of 4 (red). The results show 

23 that the cooling rate in W m' is essentially unchanged. See text for further explanation. 
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